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I. Background 


The U.S. Army Ballistic Research Laboratory (BRL) has been conducting research to 
support the design of a facility intended to subject full scale military equipment to the blast 
and thermal environment produced by tactical nuclear weapons. The design of this facility, 
termed the Large Blast/Thermal Simulator (LB/TS) is being directed by the Army Corps 
of Engineers through the Defense Nuclear Agency. An illustration of the proposed LB/TS 
can be seen in Figure 1. 



Figure 1. The Proposed Large Blast/Thermal Simulator 


As currently conceived the LB/TS, when completed, will be the world’s largest shock 
tube. A shock tube is a device in which compressed gas is quickly released from a reservoir to 
form a shock wave which travels down its expansion section. By adjusting the volume of the 
gas reservoir (also known as the driver) and the initial reservoir pressure and temperature, 
a nuclear blast wave of desired shock overpressure and weapon yield can be simulated. The 
LB/TS is designed so that a wide range of nuclear blast waves can be simulated. Shock 
overpressures in the simulator can range from 2 — 35 psi and simulated weapon yield can 
range from 1 - 600 kT. 

Conventional shock tubes usually have a single driver. The LB/TS will have a total of 
nine separate driver tubes. The downstream end of each driver tube has a converging nozzle 
and a cylindrical throat section. Located in each throat is a pair of diaphragms. The flow 
in the LB/TS is initiated through the simultaneous rupturing of the diaphragms in the nine 

1 









drivers. Some of the driver tubes, with their respective nozzles and diaphragms, can be seen 
in Figure 2, whidi is a schematic diagram of the LB/TS. 



Figure 2. Schematic Diagram of the Large Blast/Thermal Simulator 

When the diaphragms are ruptured and the gas rushes out of the driver tubes, thrust is 
generated which applies a force to each driver tube. The applied force acts in the opposite 
direction of the flow. In order to design a restraint system capable of counteracting these 
forces, the maximum expected thrust history for each driver tube must be determined. 
This report documents eflforts by the Ballistic Research Laboratory to estimate the thrust 
history for each driver tube when charged with its maximum operating pressure. In all 
the calculations performed for this study, the initial driver pressure and temperature were 
15.615 MPa and 650 K, respectively, while the ambient pressure and temperature used were 
101.325 kPa and 288.15 K. 

As previously stated, the volume of the driver tubes can be adjusted to simulate a 
wide range of blast waves. As the volume of the driver tubes is decreased, the simulated 
weapon yield for a particular set of driver conditions decreases. The volume of a driver tube 
is adjusted by positioning a device known as a hydroplug in the tube. A hydroplug, shown 
in Figure 2, is a large cylindrical device that roughly resembles a bank vault door. It is 
made primarily of steel and is equipped with metallic wedges which are expanded against 
the driver tube walls to hold the hydroplug in place. When a hydroplug is locked in position, 
its elastomeric rings provide a seal between itself and the wall of the driver, preventing high 
pressure driver gas from escaping around its perimeter. One scenario which is of concern 
to the design of the LB/TS is that in which a hydroplug should become dislodged after the 
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driver tubes have been pressurized. If this were to happen, thrust would be generated on the 
driver tube which would try to force it into the expansion section of the LB/TS. This thrust 
acts in the opposite direction of the case in which the LB/TS is operated normally. So that 
the loads associated with this type of scenario may be incorporated into the design of the 
LB/TS, an effort was taken to estimate the thrust histories in the event that a hydroplug 
should fail. 

The calculations described in this report were performed using Version 122 of the HULL 
hydrocode 2. All calculations employed two-dimensional grids with axial symmetry. The 
HULL code is an inviscid code. Therefore, forces on the driver walb resulting from fluid 
friction were not considered. All thrust histories were calculated by summing the hydrostatic 
pressure on the surfaces of the driver tubes for a given point in time. 

The gridding scheme in the HULL code uses rectangular cells. The accuracy at which 
lines that do not run horizontally or vertically are modeled is a function of the width and 
height of the cells. If the grid has very small cells, then a diagonal line, for instance, will 
be modeled with greater accuracy than it would with larger cells. The same is true for 
curved surfaces. However, the price one pays for small cells is large. In a two-dimensional 
explicit time-stepping hydrocode, the change in run time caused by a change in cell size 
is proportional to the cube of the ratio of the old and new cell widths. For example, a 
calculation with square cells 1 cm wide will consume eight times as much processor time as 
an equivalent 2 cm cell run. With this in mind, an acceptable combination of model accuracy 
and computational economy was chosen for the calculations discussed in this report. 


II. Driver Tube Layout in the LB/TS 

There are a total of nine driver tubes in the current design of the LB/TS. Four different 
driver lengths are divided among these nine drivers Two driver tubes are 14.000 m long, 
three are 18.666 m, two are 27.200 m and two are 36.159 m long. All the driver tubes 
are equipped with hydroplugs except for the two 14.000 m tubes which have hemispherical 
heads at their upstream ends. The lengths listed here represent only the cylindrical driver 
tube from the downstream end of the hydroplug to the beginning of the converging nozzle 
for those drivers with hydroplugs. In the case of the two 14.000 m drivers, the length is 
defined by the distance from the radial center of the hemispherical head to the beginning of 
the converging nozzle. 

The expansion section of the LB/TS is semicircular. The driver tubes are arranged 
so that the flow from the drivers will fill the expansion section as evenly as possible. Any 
two driver tubes of equal length are located on opposite sides of the expansion section, 
equidistant from the centerline of the expansion section. The set of three 18.666 m tubes 
has one of its tubes on each side of the expansion section with the third tube in the center 
of the expansion section. With this arrangement, no torsional moments are applied to the 
driver tube reaction pier during normal operation. The arrangement and numbering of the 
driver tub« is illustrated in Figure 3, 
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Figure 3. Location of Driver Tubes in the LB/TS, Facing the Drivers from 
a Position Inside the Expansion Section 
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III. Normal Operation 


Because there are only four unique driver tube lengths, only four thrust histories need 
to be determined. In an effort to keep the calculations as simple and economical as possible, 
only one driver tube was modeled in any one run. The computational model of one of the 
driver tubes can be seen in Figure 4. 

SETUP FOR DRIVER THRUST CALCULATIONS 
NORMAL OPERATION 

14.000 METER DRIVER TUBE MTH HEMISPHERICAL HEAD 

OnvCR SECTKM 

MSOE aAMCieX - t.a2S m CXPANSMN sechch 

LCNCIH . U.000 m W9P 6 CIMC TIW > i.S3 m 

MnAL DRIVER PRESSURE • IS6tS kPo AUBENT PRESSURE > 101J2S hPa 

M1UL DRIVER IEMPERATURE > 630 K AUBENT 1EUPEIU1UNE - 288.13 K 



Figure 4. Computational Model for Normal Operation Calculations 

The ratio of the area of the expansion section to the combined area of the nine throat 
sections is 28:1. If it is assumed that the flow from each driver tube equally fills one ninth 
of the expansion section, then the area ratio seen by the flow in going from one throat into 
the expansion section will be 28:1. Since only one driver tube would be modeled at a time, 
the area of the expansion section in the calculations is set at 28 times the area of one throat 
section. Also, in the computational model, the expansion section was made very long so that 
disturbances generated when the shock reached the downstream boundary would not effect 
the flow in the driver during the computation. 

The resulting thrust histories can be seen in Figures 5 through 8. For the most part, 
one can see that the amplitude of the pressure waves is about the same for each of the four 
thrust histories. The timing of the passing of the waves in the drivers is the cause of the 
differences between the histories. 

Each of the thrust histories starts at time equal to zero with an initial thrust of about 
10.19 million Newtons. When the flow is initiated, the gas begins to travel in the downstream 
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direction. At the same time, a rarefaction wave begins to move in the upstream direction 
into the driver. As this rarefaction wave passes over the converging nozzle, the pressure on 
the surface of the converging nozzle decreases, resulting in an increase in net thrust. This is 
illustrated by the first peak in the thrust histories at about 16 million Newtons. This peak 
is held roughly constant until the rarefaction wave meets the upstream end of the driver 
tube. When this occurs, the pressure on the upstream wall of the tube drops resulting in a 
substantial drop in thrust. After the rarefaction wave meets the upstream end of the tube 
it is reflected and begins to travel downstream where it is eventually reflected again at the 
throat and the process starts over again. It is this repeated motion of the rarefaction wave in 
the driver which gives the thrust history its cyclic pattern. Each passage of this rarefaction 
wave causes the pressure in the driver to drop. This decrease in driver pressure causes the 
underlying decay seen in the oscillating thrust histories. 

The thrust histories illustrated in Figures 5 to 8 are for single driver tubes. To obtain 
the total thrust from all nine driver tubes that would be seen by the reaction pier, the 
individual thrust histories were multiplied by the number of driver tubes for each length and 
then summed. The resulting total thrust history can be seen in Figure 9. As one would 
expect, the peak thrust is nine times the peak for one driver tube. The total thrust history 
differs from the individual thrust histories in that its cyclic pattern is less apparent. This 
is caused by the difference in frequencies between the thrust histories for driver tubes of 
different lengths. When summing the histories of varying frequencies, the decaying nature 
of the thrust history becomes more obvious. 

These thrust histories are the expected result of operating the LB/TS at its maximum 
design limits of driver pressure, driver temperature and driver volume. Lower driver pres¬ 
sures would result in lower peak thrust values. A smaller driver volume would dictate a 
quicker decay in the thrust histories. By using these maximum expected thrust histories, 
the LB/TS driver support system can be designed to withstand the expected loads under 
normal, repeated operation. 
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Thrust History for Full Seal* LB/TS Drivor Tubos 
Maximum Driver Pressure and Temperature : Normal Operation 



Time (s) 


Thrust History for 14.000 m Driver Tube; Normal Operation 


Thrust History for Full Scale LB/TS Driver Tubes 


Maximum Driver Pressure and Temperature : Normal Operation 



Figure 6. Thrust History for 18.666 m Driver Tube; Normal Operation 
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Thrust History for Full Scale LB/TS Driver Tubes 
Maximum Driver Pressure and Temperature : Normal Operation 



Time (s) 
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Thrust History for 27.200 m Driver Tube; Normal Operation 


Thrust History for Full Scale LB/TS Driver Tubes 
Maximum Driver Pressure and Temperature : Normal Operation 



Figure 8. Thrust History for 36.159 m Driver Tube; Normal Operation 




Total Thrust for Nino LB^S Drivor Tubos 
Maximum Drivor Prossuro and Tomporaturo : Normal Oporatlon 



Figure 9. Total Thrust History All Driver Tubes; Normal Operation 

IV. Hydroplug Failure 

A hydroplug is a device which is to be used to adjust the volume in the driver tubes. 
Hydroplugs are currently being used in the oil industry to block a portion of pipeline while 
repairs or modifications are being made to the pipe. However, the hydroplugs used in the 
oil industry are not as large as that required for the LB/TS. Also, the temperature and 
pressure of the driver gas in the LB/TS will be much higher than that seen in the oil 
industry application. For these reasons, a decision was made to examine the thrust load 
which would be applied to the driver tubes in the event that a hydroplug should fail after 
the driver has been pressurized. 

As previously stated, the HULL code was used to perform all the calculations in the 
study. Since HULL is a code for solving problems in hydrodynamics, it is not well suited for 
solving problems of multi-body dynamics. However, the multi-body dynamics codes that are 
available l 2 w:k sufficient hydrodynamic algorithms to provide accurate aerodynamic loading 
data. Since the fluid forces on the driver were the desired result, the HULL code was chosen 
to solve the problem. Had we been interested in the interaction between the driver walls 
and the hydroplug, a multi-body dynamics code might have been a better choice. 

The calculations of driver thrust generated from hydroplug failure were performed at 
the maximum expected driver temperature and pressure. For all of these calculations, the 
hydroplug was initially positioned at the upstream end of the driver tube. It would have 
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been possible to perform more calculations with the hydroplug in different initial positions, 
but a short schedule and the cost of computer time prohibited this. 

Of the four unique driver tube lengths in the LB/TS, only the three longest will have 
hydroplugs. The two 14.000 m tubes will be equipped with hemispherical ends. Therefore, 
thrust calculations for the hydroplug failure case were only performed for the 18.666 m, 
27.200 m and 36.159 m driver tubes. 

Three different problem geometries were modeled to achieve the thrust histories for each 
of the three different driver tube lengths. The three problem geometries differed in the way 
the hydroplug was modeled. As the accuracy of the hydroplug model improved, the level of 
complexity of the calculation increased. The results of these three geometries will now be 
discussed in detail. 


1. Simple Discontinuity 

As a first order solution to the thrust calculation problem, the hydroplug itself was not 
modeled. Rather, the driver gas and ambient air were initially in contact with each other 
at a temperature/pressure discontinuity. In this setup, which can be seen in Figure 10, no 
attempt was made to model the dynamics of the interaction between the fluids and the 
hydroplug. 

The resulting thrust histories for these three calculations can be seen in Figure 11. In 
these calculations, the thrust is initially constant at about 40.75 million Newtons. This 
thrust value is simply the net horizontal force on the driver determined by multiplying the 
initial driver overpressure by the inside diameter of the driver tube. 

As in the normal operation case, the pressure discontinuity causes the flow to travel out 
of the driver tube in the direction of the lower pressure. At the same time, a rarefaction wave 
is generated and travels in the opposite direction into the driver tube. When the rarefaction 
wave meets the nozzle end of the driver, the pressure decreases and thus the thrust decreases. 
The rarefaction wave is then reflected at the diaphragm and travels back toward the position 
of the original discontinuity. As this reflected wave passes over the nozzle, the pressure on the 
nozzle is again reduced and therefore the thrust is decreased further. Since no hydroplug was 
modeled in these calculations, no repeated passage of waves is seen in the thrust histories. 
The entire driver pressure is relieved with one passage of the initial rarefaction wave. 
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SETUP FOR DRIVER THRUST CALCULATIONS 
KYDROPLUG FAILURE CASE 

HYOROPLUG NOT MOOELED; SIMPLE CASEOUS DISCONTINUITY 



NOT TO SCALE 


Figure 10. Computational Model for Simple Discontinuity Setup 

Thrust History for Full Seals LB/TS Orivsr Tubss 
Maximum Orivsr Prsssurs and Tsmpsroturs : Hydroplug Failurs 



Figure 11. Thrust Histories for Simple Discontinuity Setup 
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2. Dense Air Mass 


The hydroplug failure calculations using a simple discontinuity which were discussed 
previously were performed as a first order estimate. These calculations were not difficult 
to perform and ran in a reasonable amount of computer time. However, it was known that 
the effect of the hydroplug on the flow was important. Therefore, further investigation was 
required. 

Due to the time stepping nature of this type of code, attempts to model the hydroplug 
as a homogeneous solid were initially avoided. The time step in a hydrocode such as HULL 
is determined by the spacing of the computational grid and the speed of sound plus the 
particle velocity of the materials in the grid. The speed of sound in a solid material such as 
steel is much greater than in air, which would result in a dramatically reduced time step. 
Therefore, this type of calculation was considered to be too time consuming to attempt. 

As a compromise, a calculation was performed in which the hydroplug was modeled as 
a dense air mass in the 18.666 m driver tube. The setup for this calculation is shown in 
Figure 12. The density of the air used to model the hydroplug was defined so that the mass 
of dense air was equal to the mass of a hydroplug. 

It was estimated that while diffusion of the air mass would occur, the diffusion might be 
small enough to accurately model the presence of the hydroplug. The diffusion of the air mass 
was observed by plotting contours of density for the computational domain at various times 
during the calculation. These density contours can be seen in Appendix A. The figures in this 
appendix show that the air mass, which served as the model of the hydroplug, experienced 
a great deal of diffusion. 

The thrust history for this calculation is in Figure 13. In this figure the thrust history 
is compared to the one obtained from the simple discontinuity case for the same 18.666 m 
driver. Figure 13 shows that the presence of the dense air mass had a significant effect on the 
thrust history. It is believed that the thrust history for the dense air mass case is a better 
physical representation of the event when compared to the simple discontinuity case. 

Like the simple discontinuity case, the thrust history for the dense air mass setup begins 
at a constant value of 40.75 million Newtons. The pressure differential across the dense air 
mass causes it to move out of the driver, in the direction of the lower pressure. When the 
air mass begins to move, a rarefaction wave is generated on its upstream face which travels 
in the upstream direction toward the diaphragm. This rarefaction wave is very similax to 
that generated in the simple discontinuity case, and therefore the initial decay of the thrust 
history is very similar to that in the simple discontinuity case. This decay is caused by the 
passage of the rarefaction wave over the nozzle. 

Because the presence of the hydroplug is modeled, the thrust history does not decay as 
rapidly as that in the simple discontinuity case. This occurs because the hydroplug does not 
allow the gas to immediately escape from the driver. After the initial decay caused by the 
first rarefaction wave, the thrust history levels out at a value of about 30 million Newtons. 
When the upstream face of the hydroplug exits the driver tube, the driver gas is allowed to 
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escape from the driver which causes another rarefaction wave to be generated. This wave 
travels into the driver tube and decreases the pressure on the walls of the driver resulting in 
a decrease in thrust. With the hydroplug now out of the way, the driver gas is allowed to 
escape with one passage of this secondary rarefaction wave which results in the exponential 
decay of the thrust history. 
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SETUP FOR DRIVER THRUST CALCUIATIONS 
HYDROPLUG FAILURE CASE 

HYDROPLUG MODELED AS DENSE AIR MASS 



NOT TO SCALE 


Figure 12. Computational Model for Dense Air Mass Setup 


Thrust Histories for 18.666 Meter Driver Tubes 
Maximum Driver Pressure and Temperature ; Hydropiug Failure 
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3. Steel Plug with Modified Properties 


In an effort to improve the accuracy of the thrust histories for the hydroplug failure 
case, a final set of calculations was performed in which the hydroplug was modeled as a 
homogeneous steel plug with modified properties. The computational model for this calcu¬ 
lation can be seen in Figure 14. As stated earlier, this approach was initially avoided due to 
the adverse effect the steel would have on the time step of the calculation. However, it was 
found that by altering some of the properties of the steel, the calculation could be performed 
without the penalty of a dramatically reduced time step. 

The density of the steel in the computational model was modified so that the mass of 
the modeled hydroplug was equivalent to the mass of a plug proposed for the LB/TS. The 
time step problem was resolved by reducing the ambient sound speed in steel in the model. 
With these modifications, improved thrust histories could be obtained within a reasonable 
amount of computation time. 

This practice of changing the properties of steel did not go without penalty though. 
The constitutive relations in this code define the sound speed in steel as the square root of 
the fraction of the Young’s Modulus divided by the density of the homogeneous material as 
can be seen in Equation 1. 



Solving Equation 1 for Young’s Modulus yields Equation 2. Prom Equation 2 it can 
be seen that by reducing the sound speed and density of steel, the Young’s Modulus is 
significantly reduced. In the calculation, the sound speed in steel was set at 500 m/s and 
the density was 4247 kgfm^. The resulting Young’s Modulus, then, is 1,062 MPa which 
is much less than the expected value of 167,042 MPa. As a result of this reduced Young’s 
Modulus, one would expect to see a great de€d of bydroplug deformation in the calculation. 
This deformation, however, was considered acceptable as long as the plug remaned roughly 
intact during the period of interest. 


£• = 0 ',, ( 2 ) 

The value of 500 m/s for the ambient sound speed of steel was not chosen arbitrarily. 
So that the steel would not control the time step in the computation, a value had to be used 
which was less than or equal to the sound speed in the surrounding air. Using Equation 3, 
the initial sound speed in the driver gas was determined to be 511 m/s. 

(3) 

In Equation 3, 7 is the ratio of specific heats of the ideal gas and is equal to 1.400, R is 
the gas constant which is 287.04 J/kg/K and T is the temperature of the driver gas which 
is 650 K. By setting the ambient sound speed in steel at 500 m/s, the computation will 
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initially be controlled by the driver gas. The steel will not control the time step until late 
in the computation when the sound speed of the driver gas has dropped due to expansion. 

As in the dense air mass case, the displacement and shape of the hydroplug were ob¬ 
served by plotting contours of density for the computational domain at various times during 
the calculation. Density contours for the 18.666 m driver tube can be seen in Appendix B. 
These figures illustrate that while the homogeneous steel plug did undergo some diffusion, it 
remained roughly intact during the period of interest. For this reason the fiow of the driver 
gas around the plug for this calculations is considered to be more accurate than the previous 
two models. Therefore, the thrust histories for the calculations which used the homogeneous 
steel plug model are assumed to be the most accurate. These thrust histories can be seen in 
Figure 15. 

The thrust histories in Figure 15 are very similar to the thrust history for the dense air 
mass case. In each, the decay of the thrust is caused by two rarefaction waves. The first 
wave originates at the upstream face of the hydroplug. The second wave is caused by the 
exiting of the hydroplug from the driver tube and subsequent acceleration of the driver gas. 
The primary difference between the thrust histories in the two setups is the rate of decay of 
the thrust caused by the first rarefaction wave. 

Since the hydroplug was modeled as a compressible substance (air) in the dense air 
mass case, the high pressure driver gas was able to compress the hydroplug and cause it to 
diffuse. This resulted in a strong initial rarefaction wave. 

In the case where the hydroplug was modeled as a solid steel plug with modified prop- 
. erties, the plug was not compressed by the driver gas and also did not experience as much 
diffusion as in the dense air mass case. For this reason, the initial rarefaction from the 
upstream face of the hydroplug did not cause as much of a decay in the thrust. The effect 
of the second rarefaction wave is nearly identical for the two cases. 

The limited diffusion of the steel plug and the fact that the hydroplug was modeled as 
an incompressible medium leads us to believe that the thrust histories for this setup were 
the most accurate of the three cases used to model hydroplug failure. For this reason, this 
setup was used to generate thrust histories for each of the three driver tube lengths which 
are intended to receive hydroplugs. 
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Figure 14. Computational Model for Steel Plug Setup 
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Maximum Driver Pressure and Temperature : Hydroplug Follure 
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Figure 15. Thrust Histories for Steel Plug Setup 
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V. Summary 


Calculations have been performed to determine the thrust resulting from the release of 
high pressure gas from various lengths of driver tubes in the proposed Large Blast/Thermal 
Simulator. In 2dl calculations, the maximum allowable driver gas pressure and temperature 
were used so that the maximum expected thrust histories would be obtained. The thrust 
histories were first calculated for the normal operation of the shock tube using four unique 
lengths of driver tube. 

Thrust histories were also calculated for the case in which a driver volume adjusting 
device (hydroplug) failed. This event would cause the thrust to act in the opposite direction 
as the normal operation case. The calculations for the plug failure case were performed for 
three different driver lengths. Several iterations were performed on these calculations, each 
successive calculation more accurate than the last. 

Using the results of this study, the support structure for the driver tubes of the LB/TS 
can be designed in such a way that the driver tubes will be sufficiently restrained during the 
operation of the facility. 
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APPENDIX A: Density Contours for Dense Air Mass Setup 
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This appendix contains the density contour plots from the hydroplug failure calculation 
which used a dense air mass to model the hydroplug. In the discussion of this calculation, 
the dense air mass was described as having experienced a great deal of diffusion as a result of 
the pressure imbalance across the hyrdoplug. The following figures illustrate this diffusion. 
In all the figures, the axis of symmetry of the driver tube runs vertically with the high 
pressure driver gas below the hydroplug and the ambient atmosphere above the hydroplug. 
The driver tubes in the LB/TS are very long relative to their diameter. For this reason the 
radial and axial lengths are drawn to different scales. 

Figure A-1 shows the setup at the beginning of the calculation. In this figure, the 
hydroplug is represented by the horizontal lines running across the page. These horizontal 
lines represent many different contours of density which are lying on top of one another due 
to the abrupt changes in density caused by the transitions from driver gas to hydroplug to 
ambient atmosphere. The following figures show the hydroplug as it is accelerated out of the 
driver tube and into the surrounding atmosphere. As the figures progress, one can see that 
as the hydroplug is accelerated it is also beginning to experience diffusion. This diffusion 
advances as in each successive picture, the contours grow farther apart. By 120 ms the 
diffusion has progressed to the point where only the bottom two contours are visible. The 
remainder of the hydroplug has spread to the point that it is no longer in the picture. 

These figures illustrate that modeling the hydroplug with a dense air mass provides 
an improvement in the calculation of thrust histories over the simple discontinuity case. 
However, the diffusion of the dense air mass was enough that further investigation should be 
performed to provide a model which more accurately describes the response of the hydroplug. 
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Figure A-1: Dense Air Mass at 0.0 ms 
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Figure A-3: Dense Air Mass at 40.0 ms 
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Figure A-5: Dense Air Mass at 80.0 ms 
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Figure A-7: Dense Air Mass at 120.0 ms 
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Figure A-8: Dense Air Mass at 140.0 ms 
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Figure A-10: Dense Air Mass at 180.0 ms 
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APPENDIX B: Density Contours for Steel Plug Setup 
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This appendix contains the density contour plots from the hydroplug failure calculation 
which used a block of homogeneous steel with modified properties to model the hydroplug. 
This calculation was performed as an attempt to improve on the results obtained from the 
calculation which modeled the hydroplug as a dense mass of air. 

As seen in Figure A-1, the system at the beginning of the calculation can be seen in 
Figure B-1. In this figure the hydroplug is represented by the horizontal lines running across 
the plot. The subsequent figures show the displacement of the hydroplug caused by the 
high pressure driver gas. By comparing plots in Appendix B with plots in Appendix A for 
the same points in time, one can see that the hydroplug, when modeled as a steel plug, 
experienced much less diffusion than the dense air mass case. The improved accuracy of the 
steel plug calculation is most evident when comparing Figure A-8 to Figure B-8. Both these 
figures represent the system at a point in time 140 ms into the calculation. In Figure A- 
8, most of the hydroplug has exited the picture and only the bottom remains visible. In 
Figure B-8, the hydroplug is not only visible, but it is still roughly intact although it has 
undergone some deformation. The limited diffusion of the hydroplug is evident in all the 
figures in Appendix B. Based on the results of this calculation it was concluded that these 
results provided an accurate representation of the loads placed on the driver tubes of the 
LB/TS should a hydroplug failure occur. 
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Figure B-2: Steel Plug at 20.0 ms 


38 











Radius cm 

60.000 ms Cycle 
18.666 METER DRIVER 


Radius cm 12 OCT 90 

1555 Problem 18.6660 
FAIL WITH STEEL PLUG 


Figure B-4: Steel Plug at 60.0 ms 
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Figure B-5: Steel Plug at 80.0 ms 
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Figure B-6: Steel Plug at 100.0 ms 
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Figure B-7: Steel Plug at 120.0 ms 
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Figure B-8: Steel Plug at 140.0 ms 
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Figure B-10: Steel Plug at 180.0 ms 
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Figure B-11: Steel Plug at 200.0 ms 
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